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ABSTRACT 


A discussion of digitizing noise and its effect on spectral sig- 
nal to noise ratios subject to the limitations of the telemetry sys- 
tem forms the main body of this study. Simulated and actual test 
data are compared with the theoretical calculations and an optimum 
working code is derived. The code used in the Nimbus B Infrared 
Interferometer Spectrometer is found to be very close to the opti- 
mum condition for the number of bits used per data word. 
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THE EFFECT OF DIGITIZING NOISE 
ON SPECTRAL SIGNAL TO NOISE RATIOS 
IN FOURIER SPECTROSCOPY 


I. INTRODUCTION 

The purpose of this study is to determine the number of bits necessary to 
represent each data point of the record sent by the Infrared Interferometer 
Spectrometer (IRIS) instrumentation such that the signal to noise ratio in the 
spectrum due to digitizing effects is kept to a minimum, and that the number of 
bits per word is such that the total number transmitted stays within existing 
telemetry limits. The study will be conducted in two parts; a theoretical descrip- 
tion of the effect, and computer simulation for a test spectrum which consists 
of emission and absorption lines. A study of thermal or random noise effects 
is also included in combination with the digitizing noise. 

The primary instrument in the IRIS experiment is a Michelson interferometer 
which produces a record from which it is possible to obtain a spectrum of a 
source by calculating the Fourier transform of the record which is recorded as 
a function of time or distance. This method of obtaining a spectrum known as 
Fourier spectroscopy and techniques for this are discussed extensively in the 
literature.®® In this study, we are primarily concerned with the types of noise 
on this record, and how the noise transforms in the spectrum. 

The two types of noise that are most common on the recording (henceforth 
called the interferogram) are (1) detector noise due to electron movement within 
the detector and which is a function of temperature, but not a function of incident 
radiation, and (2) digitizing or quantum noise which is caused by sampling the 
amplitudes at discrete levels. The latter is a function of the analog to digital 
system used, and on ground is not a serious problem since five or six place 
decimal A to D converters are available, however, on satellites a restriction on 
the accuracy is imposed by the telemetry system. For instance, with the present 
Nimbus B system, the IRIS telemetry has a transmission rate of 3.75 kilobits 
per second. For an 11 second scan and 3408 words, this yields 12 binary bits 
per word of which 4 bits used are for housekeeping. Hence, one is limited to 
sampling 256 descrete levels of voltage. The round off error due to this arrange- 
ment may be considerable even if one has a noiseless detector. This error es- 
sentially introduces noise on the interferogram of a type similar to that of the 
detector which is assumed to be white noise. 
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II. NOISE THEORY 


We now consider the question of the computation of the signal to noise ratio 
in the spectrum for a given signal to noise ratio in the interferogram. The ex- 
pression to be derived makes the assumption that the noise in the interferogram 
is uncorrelated to the signal. 

The autocorrelation theorem ® for aperiodic functions is 


I CO fOO 

f(x)f(x+r)dx = I I F(cr) | 2 da (1) 

-00 *'-00 


where F(cr)is the Fourier transform of f(t) and is related by the transform 
pairs 


F(ct) 



f ( x ) e - i2mTX 


dx 


(2a) 


and 


f(x) 



e i 2 7TCTX 


dcr . 


(2b) 


In the above notation f(x) corresponds to the interferogram or noise record 
or a linear combination of the two, and F( cr) is the complex spectrum assuming 
f ( x ) is real. From Equation (2b) the following property may be deduced: 


f(0) 



F(o-) dcr 


( 3 ) 
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which means that the amplitude of the interferogram at zero path -difference is 
equal to the total spectral intensity. Also, from Equation (1) at r = 0 we obtain 



(x) dx 



! F(o-) | 2 d a 


(4a) 


which states that the square of the interferogram equals the square of the 
spectrum, that is, the energy densities of the two are equal. Assuming finite 
limits and discrete sampling with distance or time at about the Nyquist rate 
of two samples per highest frequency (A), we can write 



(x) dx 
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f 2 (iA) 




2 


where dx-A, X = NA, cr m = l/A, dcr- i/nA, and Equation (4a) becomes 



(4b) 


The sums over i and j are the mean square values of the interferogram 
and spectrum respectively. Considering only a noise record and using n ^ and 
n g 2 for the mean square values, we can write 


n 


s 



(5) 
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where n x and n s are the root mean square values of the recorded noise and its 
transform respectively. 

We can define the signal to noise ratio in the interferogram as the peak 
signal (which occurs at zero retardation) divided by the rms noise in the inter- 
ferogram. For the spectral domain we can consider an average signal S which 
is the total spectral intensity divided by the spectral bandwidth. Referring to 
Equation (3) 


S 


r 


F(ct) der 


2o 


f(0) 

2m s NA 


(6) 


where m s is the number of resolved spectral elements in the spectral bandwidth, 
hence the signal to noise ratio in the spectrum becomes 



NAf(O) fe 
2M n n I V NX 




( 7 ) 


Equation (7) assumes that the noise and spectral bandwidths are equivalent. If 
the noise bandwidth is greater than the spectral bandwidth. Equation (7) becomes 



( 8 ) 


where M is the number of resolvable elements in the noise bandwidth. The 

n 

above assumes unity gain in the electronics hence, if we have a gain G, the peak 
signal becomes 


M°) 


f (0) G 


2m SG 
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The mean square noise on the interferogram may be considered to be the 
sum of the mean square values of the digitizing and detector noise respectively, 
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since they are uncorrelated. As the gain affects only the detector noise we obtain 


n 


i 



( 10 ) 


where n t and n 2 are the digitizing and detector noise respectively. With this 
assumption Equation (5) becomes 


n 


s 



frnr 

VNI CT m 

r n m 


(ii) 


If one observes a source with very few lines, G will be large, hence 

n s % Gn 2 /NX/M n ~ m . Using Equation (6b), the average spectral signal will be 


and 


SG 


f G (0) NA 
2M 

n 



(12) 


If is larger than Gn 2 as may be the case for a many lined or broad spectral 
source, n s - rij y<r M N/M 2 X (the case for no gain) and 



For the case where n x % Gn 2 we have 



(13) 


( 14 ) 
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In practice we attempt to keep n x less than Gn 2 . As can be seen from 
Equation (10) , a factor of two yields about a 10% error in the estimated spectral 
signal to noise ratio. 


HI. A MODEL FOR DIGITIZING NOISE 

In practice, the number of binary bits per data word limits the number of 
quantizing levels. Let K = 2 k be the number of levels and 2f ( 0 ) be the peak 
voltage. This latter is due to the fact that f(x) may be positive or negative 
with a maximum value of | f ( 0 ) I . Let V (jE 0 ) be the assigned value for f(x) if 
jE 0 < f(x) < (j + 1) E 0 j = 0, 1, ...» K - 1, where 2E Q = 2f(0)/K . This is 
shown in Figure la. 

For the sampling interval A, f(x) may in general maybe approximated by 
a straight line of slope m®as shown in Figure lb. Thus the error due to 
quantizing may be represented by 


, , . E o E o 

e(x) = mx where < x < — 

v ' m — m 


(15) 


The mean square error is then 


_m_ f E ° /m 

2E o J. E ‘ 

L 0/m 


E o 2 


( x ) dx - 3 


(16) 


for each data point. The rms error voltage is then f(0)/K V3 , hence, the signal 
to noise ratio in the interferogram for digitizing noise alone is 



1.732 x 2 k . 


(17) 


Another scheme for quantizing the interferogram is to quantize points less 
than a certain percentage of full scale to 2 k and use a coarser quantization for 
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the remaining points. For instance. Nimbus B quantizes the points of the inter- 
ferogram that are less than 10% of the maximum to 256 levels. The remaining 
points are then divided by 10 and quantized, the operation being recorded by one 
bit of the telemetry data word. The mean square noise due to this is 


y 


(V ) 2 e » 2 

3 + c 1 ~y) t 


where y is the percentage of total points in the interferogram that are digitized 
with 2E 0 ' as the step size. The coarse quantization step size is 2E 0 . Since 


f(0) 

2 k 


and 


pf(Q) 

2 k 


P£ 1 


(18) 


the rms noise becomes 



"I) + 1 


and the signal to noise ratio due to digitizing is 


1.732 • 2 k 

yV (p 2 ~ l) + 1 



(19) 


where the denominator is referred to as the increase factor. Since y and p are 
always less than one, an increase in S/N ratio results with the same number of 
bits. As an example if 90% of the words are less than 10% of the maximum value 
we have an increase of a factor of 3 over equation (12). Figure 2 shows a plot 
of the increase factor as a function of p. 
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IV. TEST RESULTS 


A model spectrum® consisting of emission and absorption lines was used 
to generate a symmetric interferogram which was then quantized to different 
levels and retransformed. Figure 3a shows the correct spectrum; Figures 3b-3d 
show the spectrum of interferograms normalized to 256, 1024, and 2048 levels 
(k = 8, 10, 11); Figure 3e shows the spectrum computed by applying the Nimbus B 
scheme. The noise in this spectrum lies somewhere between that of Figure 3c 
and 3d. A look at the interferogram revealed that y£.98. From Figure 2 the 
(S/N) increase is about 5 and corresponds to 2500 sampled levels. 

Figure 4 shows the total percentage of points falling below a given per- 
centage of the maximum voltage accepted by the IRIS system. The data is the 
average of 13 interferograms looking at the atmosphere obtained by the University 
of Michigan balloon experiment. An interesting feature of this is that the 
Nimbus B acquisition system is nearly optimal as about 98% of the data points 
are less than 10% of the maximum. 


V. SUMMARY 

The formula for determining the signal to noise ratio in the spectral domain 
is found as a function of the signal to noise ratio of the interferogram and is 
given by equation 8. A general technique for comparing various digitizing 
schemes is given, and it is found that the Nimbus B IRIS scheme is near the 
optimum. In general, with proper digitizing codes, quantization noise is not a 
serious consideration in the signal to noise ratio of the spectrum. 
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Figure 3-Test spectrum with inteferogram quantizing effects. 
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Figure 3-Test spectrum with inteferogram quantizing effects 







